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A great deal of research has recently been conducted to solve the subsyn- 
chronous rotor vibration problems In high-performance turbomachinery. Partic- 
ularly, the destabilizing effect of the labyrinth seal on compressors or 
turbines has been Investigated for many years (refs. 1 to 9). In spite of many 
efforts the dynamic effect of the labyrinth seal had not been fully determined 
from qualitative and quantitative points of view. But from our theoretical and 
experimental work, we have determined completely the dynamic characteristics of 
the labyrinth seal. 

This paper presents the results of recent theoretical and experimental 
works. 

We developed a theoretical study and a numerical calculation program to 
obtain the dynamic coefficients based on Iwatsubo's perturbation method 
(ref. 3) and Jenny's tangential momentum effect evaluation method (ref. 9). 

The simplified formulation was programmed for practical design use. Qualita- 
tive and quantitative evaluations of the computer program have been done In 
several published works. Our experimental study also evaluated damping coef- 
ficients and considered Inlet swirl effects. 

Experimental studies on the labyrinth seal have been performed to Improve 
blading efficiency and to analyze rotor dynamics. For example, the basic lab- 
yrinth seal test was done In 1970 to verify Alford's theory, and static and 
semistatic tests were performed to Improve design, to reduce leakage, and to 
evaluate cross-coupled stiffness. In 1984-1985, to confirm the phenomena, the 
theoretical analysis of dynamic coefficients, and the swirl effect of the lab- 
yrinth seal, we continued seal dynamic model tests. This paper presents pri- 
marily the results of the dynamic test. 


SYMBOLS 

a,b displacement 

C peripheral velocity 

f cross section of seal chamber 

g gravity acceleration 

h strip height Pft£CEDlf*G PAGE SU*KK KGT 
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L length of labyrinth seal 

1 length of strip pitch 

n ratio of specific heat 

P pressure 

q mass flow rate In axial direction 

R gas constant 

R $ radius of labyrinth seal 

T absolute temperature of gas In seal 

t time 

U' length of acting surface of shear (stator) 

U" length of acting surface of shear (rotor) 

u peripheral velocity of labyrinth seal, R s »ft 

w peripheral unit length, R s 

a radial clearance of seal 

e a ,e b angu1ar displacement 

friction coefficient (stator) 
x" friction coefficient (rotor) 

v strip flow coefficient 

p density of gas 

t 1 friction shear stress of stator surface 
t" friction shear stress of rotor surface 

ft rotating speed of rotor 

o) whirling speed of rotor 

Subscripts : 
a outlet 

e entry 

F strip 

1 seal chamber number or strip 
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x axial 

Z strip number 

* steady state 


THEORETICAL STUDY OF DESTABILIZING FORCE CAUSED BY LABYRINTH SEAL 

To Investigate the destabilizing force caused by the labyrinth seal, an 
analytical model of the labyrinth seal was established for calculating eight 
dynamic coefficients (four stiffness coefficients and four damping coeffi- 
cients) considering Inlet swirl effects. 


Modeling the Labyrinth Seal 

In the flow model of the labyrinth seal Kostyuk Introduced one peripheral 
velocity variable C In the core flow of each labyrinth chamber and developed a 
simple equation (ref. 5). The developed analytical method uses the modified 
Kostyuk equation on the labyrinth seal shown In figure 1. 

The following fundamental equations are developed for the differential 
element of unit length showed In figure 2: 


Mass Flow Rate Passing Through Strip 


q? = yf-sf-CPf-! - P?) 


( 1 ) 


Mass Flow Rate Rectified In Chamber 


2 irRsFi>qi = 2 trRsi # qei 
2 irRsFi+l.qi+i = 2 irRsi.q a i 


( 2 ) 


Continuous Flow Rate In Chamber 


9 (Pifi) , f - 9 (PiCi) 

3t 1 awi 


"P (lai ” *Iei^ ® 


(3) 


Circumferential Momentum In Chamber 


9(Pi fiCi) , , 9( PiCi 2 ) , , „ „ x . 9Pi 

at + fi ^ + (< lai C ai- ( lei C ei) + T i U i-^iUi - _f i'— 


3wi 


9wi 


( 4 ) 


Equation of State 


Pi = gPi’Ri'Ti 
Pi = pi = Const. 
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These equations were established for each strip and chamber. And the Inlet and 
outlet conditions of the seal were given as follows: 


To — P*o — ?e P z - P*z — Pa 


To = T e 


Pz = T a 


Cq — C*o — C ( 


( 6 ) 


Method of Solution 

To solve equations (1) to (5), we applied Iwatsubo's method (ref. 3), that 
Is, the perturbation llneallzed method, as follows. The following nondlmen- 
slonal variables n, f, and ^ were Introduced as 

Pi = P*iU + £i), Ci = c*id + m) 

Cei = Ce;'<i (1 + Hei)> C a i = Ce , i+l = Ce&i+l (1 He> i+1 ) (^) 

q i = q*i*(l+ Ci), 5 i = 6* i( l + i) 


and, assuming that the rotor Is whirling along an elliptical orbit, 
Is represented as 


4>i = 



COSU)t *COS l P 


+ 



sinoit -sin 1 ? 


( 8 ) 


Rotor displacement a*,b* and angular displacement e a ,6b have 
the following relation: 


i-l 

a* + 6 a • X 2--; 

j=l J 


bi = b* - 


i-l 

Bb'S £j 

j=l J 


( 9 ) 


Then these equations were divided Into the steady-state equations and the 
dynamic-state equations shown In table I. As the number of variables was 
greater than the number of equations, the following two assumptions were made: 

(1) Steady-State Tangential Momentum Parameter K s 

The parameter K s , suggested by Jenny (ref. 9), Is defined as follows: 

Ce-,vi ” C a *i ” Ks*(C e *i ~ C*i ) (10) 

This parameter Is the one most Important to the destabilizing force and depends 
on the labyrinth seal geometry. 

(2) Dynamic-State Tangential Momentum Parameter K 0 

In the dynamic state, a parameter Kg, different from Jenny's K s 
parameter (ref. 9) as 
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Hei * K 0 .r|i 


01) 


These steady state and dynamic state differential equations yield to the alge- 
braic linear equations with eight coefficients K xx , K xy , K yx , K yyt C xx , C xy , 
Cy X , and C yy by Iwatsubo's method (ref. 3). 

And assuming that the rotor Is at the center of the labyrinth seal, the 
coefficients satisfy the next condition 


° K xx = K yy> K xy = “ K yx> c xx = c yy» c xy = ” c yx (12) 


Numerical Analysis and Comparison Between Theory and Published 

Experimental Results 

Two experiments on labyrinth seal destabilizing force have been published: 
Wright's (ref. 1), on the effect of bore taper; and Benckert's (ref. 12), which 
clarified the effect of entry swirl. First, the analytical results of using 
the preceding method are compared with Wright's experimental results. The con- 
figuration of the seal Is shown In figure 3. The calculated and measured 
dynamic coefficient data are shown in figures 4 and 5. The calculation was 
performed with respect to the experimental data on the effects of taper bore. 
The taper bore effect Is summarized In table II. 

The second step compares the calculated results with Benckert's experi- 
mental results for the full labyrinth seal. As shown In figure 6, the calcu- 
lated results and Benckert's experimental results are compared using Benckert's 
nondlmenslonal variables K* 0 , E*q as follows. 


* _ Kxy 

K o (P Z -P 0 )*R S -L (13) 

6 

P O ^ x / -T-) p i Po _ 2 v 

2 ^ ~ ”^axo^ ( 14 ) 


The calculation and experiment have a good agreement. The calculated entry 
swirl effect Is also shown In table II. 


EXPERIMENTAL STUDY OF SMALL LABYRINTH SEAL MODEL 

A small labyrinth seal model was tested to qualitatively confirm labyrinth 
seal dynamics. The experimental model Is shown In figure 7. The model casing 
had four nozzles attached to the annular chamber of the labyrinth seal In the 
tangential direction (ref. 7). The Inlet swirl could be alternated by nozzle 
selection for each test condition. The dimensions of the model labyrinth seal 
are summarized In table III. The model was designed to demonstrate the occur- 
rence of whirl at relatively low pressure. 
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The main test Items are summarized In table IV. The tests measured system 
damping for each test condition. The effect of shaft rotation Is very small 
because of the size of the model and the limit of the rotating speed. There- 
fore most tests were preformed In nonrotating conditions. System damping was 
measured by perturbing test working conditions. The free vibration decay was 
measured for each test. The following results were obtained from this series 
of tests. 


Effect of Labyrinth Seal on Rotor Stability 

Figure 8 shows the typical test results for the original straight seal. 
System damping varied according to nozzle Inlet pressure. The nozzle Inlet 
pressure represents the seal Inlet swirl velocity. The seal Inlet pressure was 
about one-half of the nozzle Inlet pressure. 

The measured damping ratio tended to Increase up to 0.2 kgf/cm 2 , to 
decrease as pressure Increased, and to fall Into the unstable region for pres- 
sure over 0.5 kgf/cm 2 . The vibration waves In figure 8 clearly show the 
change of system damping. 


Effect of the Labyrinth Seal on Damping 

Figure 9 shows test results at the no-swirl condition for the original 
straight seal. The damping Increased with Inlet pressure and the natural fre- 
quency slightly decreased. This shows that the seal has a direct effect on 
damping. 


Effect of Tapered bore 

Figure 10 shows the test results for the simplified tapered-bore seal. 
The clearances were changed for half the number of seal fins so that the seal 
would simulate both a convergent and a divergent seal. For this model the 
convergent seal showed more stable characteristics than the divergent seal. 
However, the differences between them were very small. 


Effect of Swirl Breaker 

To reduce the destabilizing effect of Inlet swirl, two types of swirl 
breaker were tested. The one had radial bypass holes and the other had antl- 
swlrl bypass holes. The results (fig. 11) show a significant Increase In 
stability limit for both cases, and with the antlswlrl breaker the stable con- 
dition could be maintained to about four times the Inlet pressure. 


Comparison of Test Results with Calculated Results 

Figure 12 shows the nondlmenslonal destabilizing effect (by Benckert's 
method) for swirl test results and analytical values for conditions associated 
with the original straight seal model and with tapered-bore seal models. The 
figure shows fairly good agreement between the theory and the experiment for 
both the qualitative and quantitative points of views. 
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EXAMPLE OF APPLICATION 


The method of calculating labyrinth seal dynamics was applied to solve 
compressor vibration problems (ref. 10). 

When the compressor was replaced by a new machine with Improved perform- 
ance the machine experienced severe unstable subsynchronous whirl over 90 per- 
cent load. The stability characteristics were analyzed by the Mitsubishi rotor 
dynamics program (ref. 11). The middle of figure 13 shows the stability graph 
of this rotor system; the graph Includes the labyrinth seal destabilizing 
effects calculated by this work. 

After lengthy discussions of the analytical results and the observed phe- 
nomena, we decided on a countermeasure, the Installation of a damper bearing. 

A one-day shutdown of the compressor allowed the damper bearing to be Installed 
without unbolting the compressor casing. When the compressor was run with the 
damper bearing, the subsynchronous vibration completely disappeared. The top 
and bottom figure of figure 13 compare vibration records from before and after 
damper bearing Installation. 


CONCLUSIONS 

Our theoretical and experimental study of the destabilizing force of the 
labyrinth seal confirmed the following dynamic characteristics: 

1. The unstable vibration phenomena of labyrinth seals are clearly 
demonstrated by a simple model rotor system. 

2. The existence of the damping effect In labyrinth seals Is confirmed In 
the absence of Inlet swirl. 

3. For this model the tapered clearance of the labyrinth seal has little 
effect on the destabilizing force. 

4. The special swirl breaker showed a reasonable reduction of the desta- 
bilizing effect of the labyrinth seal. 

5. Application of the results of the stability analysis gave a reasonable 
Interpretation for actual turbomachinery vibration phenomena. 
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TABLE I. - LINEALIZED EQUATION OF LABYRINTH SEAL 


Circumferential Momentum Equation 


-f|i i£i + f*i 9 ni. + i*L , c*i\ agi 2 &ic*i ani 

n 9t 9t Rsi V c*i n / 34) + Rsi 8<P 


_ 8 Rj -TiP li+ml-n <5 li+iCesvj+lRsFi+l 
^e*iP>viC*iR 2 si 


5i+l 


i r~ SRiTiP^i s 2 x 2 - o 9 9 o 

+ L qeydC*iRlr Ui+l 65 Vi+lCe * i+lR SFi+l+ ]X±& *iC e*iRSFi ) 


+ |{XiUi | C*i | - A iUi j ui - C*i | (~7 " 1 ) } ] Ci 

gRiTiPfti-l Vli<5ftiCe>'<iRsFi r gRiTiqe-.':iCe-.':i+l 

qe*iP*iC*iR s i * ^ P*iC*i ei+1 


+ {XiUi I c*i I + Ajul I m - c*i I >n i - 8RiTiqe>viCeAi n ei 

1 P*i C*i 


gRiTi 


2 i p * i C * iRg i ^i + l 6 *i + lCe*i+lR^Fi+l(P*i - P * i+l) - y - 6, d Ce, ; iRl Fi ( P^-l" p|i) } 


{ a*Cos ( 4> +wt )+a* cos ( ! fnit) -b*cos(’P +ojt )+b* cos ( ip -^t ) } 


cu it L 


+ 2 { " a * s l n ( ^+wt)+a*sin( l P -ojt)+b*sin( ; P + w t)+b*sin( P- W t)} 


+ 2qe*iP^iC*iR|i ^ 2i+l6 *i+lCe*i+lR| F i + l(Pli-Pli+l) ±lj 

- yi <5*iCe*iRsFi(P*i-l- P*i) S^j} 

j=l 

{0 a COs( l P+(i)t)+ 0 a cos(lp -OJt )- SbCOsC-P +<jjt)+0bCOsOP ~ U)t ) 

+ ^ j “ S i 8i ) { -0 a sin(V+cDt )+0 a sin(^-w t)+0bsinOP +u)t )+0bsin( l P -wt) } 
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TABLE I. - CONCLUDED. 


Steady State Equation 


° Mass Flow Rate Equation 


qi 2 

Ui 2 6i 


= p - p 

. 2 r . 2 L l-l r i 


o Continuous Equation 

RSFiqi = RSiqei = RSFi+lqi+1 


Circumferential Momentum Equation 


qe*i(Ce*i+l - Ce*i) + j P*i|C*i|C*i 


1 AlUi 

2 gRiTi P * i l ui ' C*i.|( u i - C*i) = 0 


Dynamic State Equation 


° Continuous Equation 


^Si , C*if fti 9£i C*if*i dr]± gRiTi P*i+lP i+l6 *i+lRSFi+l _ 

n at Rs-i-n Pci.n a ID ~ _ 2 _ 5 1+1 


9t Rsi*n Rsi*n 3ip 


gRiTi P*i ,22 2 


( y i+1 5 *i+lfc Fi+1+ U i 6 *i RS Fi ) - 


%*iR si 


qe*iP*i R si 

gRiTi P*i-iy i 5 *iRsFi _ . 


qe*iP*iR si 


u gRiTi 
2qe*iR*iR si 


{ y i +1 <$ *i+lR S F i +1 ( P *i - p Ai +1 ) - y i 6*iRSFi ( P*i- 1 - P*i ) } 


{a^cosC^+wt )+a*cos(^-out )-b*cosOj)+out)+b*cos(ip -out) } 

* ~ { “ a * s in ^ tout) +a* s in (l -out)+b*sin(ip+out )+b*sin(ip -out)} 


SRi Ti ,2 2 , 2 2 i 2 22 2 i~l 

• p ■ p 2 * ^i + l^*i+lRsFi+l(P*i " P*i+l) ~ y i & * iR SF i( P* i - 1”P * i) } 

Zqe*iP*iR si 


{BacosOP +u)t)+6acos( ip-out) -0bcosOP+oi)t)+0bcos(ip-oot) } 


+ 2 * (S£j-Si£i) {- 0 a sin( 1 P+out)+easin(^-ajt)+0bsin(ip4xDt)+0bsin(^-( J jt ) } 
j=l 
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TABLE II. - SUMMARY OF CASE STUDY 



Tapered bore seal 

(Diverging - straight - converging) 

Entry swirl 
(backward - forward) 

K xx 

► Rigid 

- ► A little rigid 

K xy 

► Destabilizing 

for forward swirl 

(-) (+) Destabilizing 

* for forward 

swirl 

c xx 

► Stabilizing 

Positive and almost 
Independent of entry swirl 

c xy 

Rigid for 

forward swirl 

Positive and almost 
Independent of entry swirl 


TABLE III. - SPECIFICATION OF TEST MODEL 


Seal diameter, mm 100 

Seal radial clearance, mm 0.25 

Height of seal fin, mm 2.75 

Pitch of seal fin, mm 4.00 

Numbers of fins 15x2 

Inlet pressure, atm, absolute 1 to 3.2 

Discharge pressure, atm, absolute 1 

Critical speed, rpm 930 


TABLE IV. - TEST ITEMS AND OBJECTIVES 


Test Items 

Objective 

Original 

Swirl effect 

Effect of Inlet swirl 

Rotation effect 

Effect of rotation of rotor 

Clearance effect 

Effect of seal clearance of same configuration 

Tapered 

clearance effect 

Effect of convergent and divergent clearance 
configuration 

With swirl 
breaker 

Swirl breaker 
effect 

Effect of specially designed swirl breaker 
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Seal Type 

C 2 /C, 

Ci (mm) 

C 2 (mm) 

Diverging 

Straight 

Converging 

1.4973 

1.0 

0.6642 

0.1311 

0.1585 

0.1915 

0.1963 

0.1585 

0.1272 


Outlet Pressure : 1.076kgf/cm 2 


Figure 3. - Seal configurations by Wright's 
test (ref. 12). 
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Figure 4. - Comparison of calculated and measured Figure 5. - Comparison of calculated and 

stiffness coefficients. measured damping coefficients. 
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Figure 6. - Comparison of calculated cross-coupling coefficient with that 
measured by Benckert's test (ref. 1). 


Figure 7. - Test model. 
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Figure 8. - Typical dynamic test results for straight seal. 
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Figure 






ANALYSED 

STABILITY MARGIN kgf/cm 


0.4u p-p 


OOHz 


lonoa /Z ))))))))) 

WITH DAMPER |p >>>))) n g 
BEARING 9 stage Compressor 

5000 ANALYSIS 

NO DAMPER BEARING 

MEASURED STABILITY LIMIT 

' 5000 10000 15000 

UNIT OUTPUT kW ^ 


Load Down 


Load Up 












Figure 13. - Example of application of damper bearing to centrifugal 
compressor. 









